Nonlinear Sigma Model Analysis of the AFM Phase Transition of the Kondo Lattice 
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We have studied the antiferromagnetic quantum phase transition of a 2D Kondo-Heisenberg square 
lattice using the non-hnear sigma model. A renormalization group analysis of the competing Kondo 
- RKKY interaction was carried out to 1-loop order in the e expansion, and a new quantum critical 
point is found, dominated by Kondo fluctuations. In addition, the spin-wave velocity scales loga- 
rithmically near the new QCP, i.e breakdown of hydrodynamic behavior. The results allow us to 
propose a new phase diagram near the AFM fixed point of this 2D Kondo lattice model. 

PACS numbers: 



The physics of heavy- fermion metals, systems that 
contain both localized /-electrons and conduction elec- 
tron bands, has been of great interest in the strongly- 
correlated community for several decades. Initial work 
by Doniach [1] pointed out the competition between the 
Kondo {Jk) and the RKKY {Jh) interaction, and the 
quantum critical point (QCP) associated with the tran- 
sition from an antiferromagnetic metal (AFM) to a para- 
magnetic (PM) state as Jk and Jh vary. Work on the 
two-impurity model [2, 3] showed that the RKKY and 
Kondo coupling are not mutually exclusive either for fer- 
romagnetic or antiferromagnetic interactions. In fact, a 
correlated Kondo effect is the rule rather than the ex- 
ception, with the impurities partially compensated by 
the conduction electrons, and partially by each other 
(for AF interactions). In 1-D, there have been various 
numerical [4] and analytical [l, 0] studies. However, 
the physics of higher dimensional systems and the corre- 
sponding quantum phase transition (QPT) is much less 
well-understood. 

A commonly assumed scenario is that the /-electrons 
delocalize and are included in the Fermi surface, resulting 
in a large Fermi surface as given by Luttinger's Theorem 
0, and the QPT is of the spin density- wave (SDW) 
type. This QPT has ty pica lly been understood in the 
Hertz-Millis approach [9|, |lO| , and the physics of delocal- 
ized /-electrons forming a heavy fermi liq uid have been 
well- understood in a large-N approach [ll[[l2|,[l3. How- 
ever, experimental studies have shown that the Hertz- 
Millis picture may be inadequate for describing the be- 
havior of the system near the QCP [3, An alterna- 
tive picture of the AFM to PM transition is that the /- 
electrons do not delocalize and become part of the Fermi 
surface, hence the Fermi surface is small. This picture 



has been proposed by several groups lly, , and stud- 



ied within the DMFT approach p^, [18| . Exotic ground 
states have also been proposed for Kondo lattice systems 
I9I . [20! . [2H that include fractionalized quasi-particles. 

We study the effect of the Kondo coupling on the anti- 
ferromagnetic phase transition, and also see if the various 
exotic states can be obtained from a more microscopic 



approach. We assume that the AFM state has a small 
Fermi surface, with no "hot-spots" spanned by the Neel 
vector (7r,7r), as shown in Fig. [TJ To study the AFM 




Fermi surface is not nested by Q 



FIG. 1: Diagram of Fermi surface with no "hot" -spots nested 
by Q = (7r,7r). The solid black lines indicate the Fermi sur- 
face, and the dashed black lines indicate the magnetic B.Z. 



QPT of the two-dimensional (2D) square Kondo lattice 
with a small Fermi surface, we map the system onto a 
non-linear sigma model (NLSM) coupled to conduction 
electrons. We then obtain an effective Lagrangian by in- 
tegrating out the fermionic degrees of freedom, and carry 
out a renormalization group (RG) analysis to one-loop 
within the e-expansion. The calculations and results are 
described in the following sections. 

As we shall show, there is a new QCP, with logarithmic 
scaling due to the Kondo interaction that plays an essen- 
tial role in the QPT. The Kondo interaction strongly af- 
fects the spin- wave fluctuations near the QCP, and leads 
to logarithmic scaling of the spin- wave velocity, implying 
localization of the spin- waves due to the Kondo effect. 
Based upon the RG results, we also propose a possible 
phase diagram near the AFM state. 

The Hamiltonian for the Kondo-Heisenberg model is 
written as: 
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where J^K > 0, are the Heisenberg and Kondo cou- 
phng respectively, and Sc{ri) = c] ^aafBCi^a- Since we are 
interested in the AFM phase, the Heisenberg term can be 
mapped onto the non-hnear sigma model in a well-known 

manner [22]. The local spin. Si = (— l)*n+ -i-., has both 

\s\i 

a Neel component, n, and a small ferromagnetic moment, 
/ oc drfixn. The conduction electrons couple to both the 
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Neel vector, Sc{ri) • n, and also the ferromagnetic part 
of the local spin, Sdu) ■ I. Upon integrating out the 
electrons, we obtain two 4th order terms, one due to the 
coupling to the Neel component and the other to the fer- 
romagnetic part. As we have assumed a Fermi surface 
that is not nested by the AF wave vector, (7r,7r), there 
is a kinetic energy gap for the former term. Similarly, 
there are no low-lying fermionic excitations that couple 
to ft to give rise to Landau damping; hence it remains a 
z = 1 theory. The effective Lagrangian in terms of the 
Goldstone modes, tt^, where a e x^y is. 
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Fig. [2] shows the Feynman diagrams for the spin- wave 
propagator and the two interaction vertices. 
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FIG. 2: Feynman rules for the spin- wave propagator and two 
vertices. The third diagram is the vertex due to the Kondo 
interaction. 



The spin-wave coupling, g 



effective Kondo 



coupling, gj. = j{ — ) N{Ef)^ spin- wave velocity, c 



^2dJSa, and spin- wave stiffness, ps = JS'^o?'~^^ are de- 
rived in terms of microscopic quantities: the Heisenberg 
coupling J, Kondo coupling, lattice spacing, a, spin 
length, iS, electronic DOS at Ep^ N{Ef), and Fermi ve- 
locity. For convenience, we have written k = (cj, k). 
Following methods similar to Brezin et al. 
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we define a renormalized theory at a momentum scale /i, 
and invariance of the theory under a change of /i gives 
the following Callan-Symanzik (C-S) equation. 
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We carried out a renormalization group calculation to 



d_ d N 



N 

T 



l-loop order within the e-expansion, and the system is 
taken to be (1 + (1 + e))-dimension. The l-loop correc- 
tions to the spin- wave propagator, G(/c), and the Kondo 
coupling vertex, F^^^ are calculated. The diagrams for the 
spin-wave propagator are shown in Fig. [3l and there are 
11 diagrams for the Kondo vertex which are not shown 
for convenience. We then define the following renormal- 
ization factors to absorb the divergences: the spin- wave 
coupling renormalization Z^, spin- wave velocity renor- 
malization Zc, and the Kondo coupling renormalization 
Zg^. A straightforward perturbative calculation shows 
that the spin-wave velocity is renormalized, as expected 
from the Kondo term which breaks Lorentz-invariance. 
The wave-function renormalization Z^^ is obtained from 
the l-loop correction to {(j{x)). We then obtain the fol- 
lowing beta functions. 
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FIG. 3: The two diagrams that contribute at one- loop order 
to the propagator. The left diagram is the usual spin- wave 
correction, and the right diagram is the correction due to the 
Kondo interaction. 



From /3g and we obtain the following RG-flow 
phase diagram, Fig. [H which shows that the unstable 
fixed point in the NLSM has shifted to a new tetra-critical 
point at ^ 9. lire and ou ^ ^ 10.27r6^. 




We study the nature of the tetra-critical point, and 
show that it is an infinite-order phase transition governed 
by a logarithmic singularity. We solved Eq. [3] and [H near 
the QGP, and obtained the correlation length, spin-wave 
velocity and propagator. The C-S equation is solved in 
a standard manner using the method of characteristics, 
and is given by 

ri'\u,c\klg,gk) = exp-^.^^^^^^'^^-^- 

Tl^\ij,Cr\k\,gr,gk,r) (5) 

where, c^, gr and gk^r are solutions of Eq. HJ and p = 
Y^cJ^ + d^l^ is the momentum scale we are interested in. 

Solving the Hessian of Pg and gives two negative 
eigenvalues, showing that the tetra-critical point is an 
unstable fixed point, and the dominant eigenvalue lies 
along the separatrix between regions I and IV. Since the 
three coupled (3 functions are difficult to solve in general 
near the critical point, we will solve it along the direc- 
tion of the dominant eigenvalue, where we know ^ as a 
function of gk and c. Linearizing the /3 functions we are 
able to obtain the following solutions. 
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FIG. 4: Plot of RG flow. 

Fig. m shows four possible phases within reach of the 
e-expansion. In region I, both g and g^ flow towards 
zero, i.e. a stable AFM phase as expected. In region II, 
g flows towards zero and gk flows away to a large value, 
indicating that the RKKY interaction between the spins 
is relevant. Since the RKKY coupling of the ferromag- 
netic moments is relevant, we expect a spiral phase with 
the ferromagnetic part of the spin having a wave- vector 
of 2/ci?. A recent spin-polarized STM study of a lattice 
of manganese atoms on a tungsten surface shows a spiral 
phase with a wavelength of = 12nm [25[. The Fermi 
wavelength is Xp = 20 nm, and 1/2 Aj? ~ Ag, indicating 
that it may be an experimental realization of the spiral 
phase of our theory (Region II). 

In region IV, g flows away to some large value, and gk 
flows to zero, which is basically identical to the quantum- 
disordered phase in the NLSM. Additional studies 2^, 27| 
have shown that proliferation of topological excitations 
drive the system into a valence-bond solid state. In re- 
gion III, both g and gk flow away to some large value, 
indicating a Kondo driven paramagnetic phase where the 
Kondo coupling is relevant and the spins are disordered. 
A possibility would of course be a heavy fermi liquid PM 
state, but more exotic Kondo-stabilized spin liquid states 
as mentioned above are also possible. 



A is the cutoff of the bare theory with the unrenormal- 
ized parameters. Here y is deflned as gk = gk,c{^ — y)^ and 
yo measures the initial distance from the critical point. 
The function ProductLogix) is also known as the Lam- 
bert function, W{x). From Eq. [6l we obtain the correla- 
tion length, ^ when y{ii) ~ 1, and ao ~ A~^. 
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This shows that we have an inflnite-order phase transi- 
tion that is dominated by a logarithmic singularity. The 
phase transition due to spin- wave fluctuations seen in the 
NLSM is now dominated instead by Kondo fluctuations 
of the spins. This clearly shows that local Kondo fluctu- 
ations play an essential role in driving the QPT, and is 
reflected in the Kondo interaction being relevant for the 
two phases found in region II and HI. 

We also solved for the spin-wave propagator near the 
QGP using Eq. \5\ which gives. 
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Note that in Eqns. [6] -[HI the exponents are numerical 
values that depend on ai and a2; the actual values are 
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h ^ i ^ 4^ ' which we closely round off in the 

exponents. 

Eq. [6] shows that for small enough momentum, p <C 
(e + 2)~57A, the spin wave velocity scales logarithmi- 
cally, which can be seen as a cross-over length-scale from 
Gaussian behavior to critical behavior of the spin waves. 
Similar behavior occurs in the Grinstein-Pelcovits renor- 
malization of the elastic constants in smectic crystals [28] . 
Therefore for a small enough region close to the QCP and 
small enough momentum, the spin-wave velocity scales 
logarithmically as, 

c(m) - Ji— (9) 

Thus for asymptotically small momentum, the spin- 
wave velocity vanishes logarithmically, indicating break- 
down of hydrodynamic behavior. This implies that the 
spin waves become increasingly localized near the QCP, 
and that localized Kondo-induced spin flips that fluctu- 
ate temporally may also be critical modes at the QCP, 
which may be a mechanism for driving the system into a 
heavy fermion or new exotic phase. Similar ideas of local 
quantum criticality have been proposed by other groups 
[lil . H^, and a similar study using Shankar's fermionic 
RG method found a Lifshitz transition [30]. 

In conclusion, we have studied the antiferromagnetic 
quantum phase transition of a Kondo lattice with a small 
Fermi surface. Since the Fermi surface has no "hot"- 
spots nested by (7r,7r), there is no Landau damping of 
the spin waves, and the theory remains z = 1. The effec- 
tive theory contains a new Kondo driven interaction that 
signiflcantly renormalizes the spin wave velocity, and is 
essential in driving the quantum phase transition. 

The RG results suggest a possible phase diagram with 
two disordered phases, one with the Kondo coupling rel- 
evant, and the other with the Kondo coupling irrelevant, 
where exotic phases with fractionalized quasi-particles 
may be found. One result of our calculations is that 
the Kondo coupling is irrelevant in the magnetic phase, 
but is relevant in one of the disordered phases, which is in 
agreement with some of the experimental results showing 
a change in Fermi surface size (l^ . 

The key result of our RG calculations is that the QCP 
of the NLSM is now shifted to a Kondo-driven inflnite- 
order QCP that is characterized by a logarithmic sin- 
gularity. The velocity of the critical modes are loga- 
rithmically slowed down due to Kondo fluctuations with 
the conduction electrons near the QCP, leading to lo- 
calization at sufliciently small momenta. This clearly 
shows that local Kondo physics plays an essential role 
in the quantum phase transition, as opposed to the stan- 
dard view of a spin wave driven QPT that is seen in the 
NLSM. The logarithmic scaling of the spin-wave veloc- 
ity is a quantum analog of a similar logarithmic scaling 
of the elastic coefficient seen in a classical Lifshitz tran- 
sition. The appearance of localized critical modes at a 



QCP, termed local quantum criticality, is a topic of great 
current interest. An outstanding issue is the effects of 
topological excitations on the QCP that we have found, 
especially in the presence of Kondo coupling to conduc- 
tion electrons. Our work could be of relevance to heavy 
fermion systems that display effective 2D-like transitions, 
and also to STM-engineered 2D Kondo lattices. 
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